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ABSTRACT 



Context. V838 Mon erupted at the beginning of 2002 becoming an extremely luminous star with L = 10' Lq. Among various 
scenarios proposed to explain the nature of the outburst, the most promising is a stellar merger event. 

Aims. We attempt to obtain information about the structure and evolution of the object in the decline from the 2002 eruption. 
Methods. The results of spectroscopic observations of the object obtained in October 2005 with the Keck/HIRES instrument, 
presented in detail in Paper I, are analysed and discussed. Modelling of the observed line profiles has been used to constrain the 
physical parameters of the system. 

Results. The kinematics of the atmosphere of V838 Mon is very complex. Our analysis of the molecular bands and the P-Cyg profiles 
of atomic lines shows that the object loses matter with a velocity of up to 215 kms"' and a rate of 10"' - 10"' Moyr"'. In the 
profiles of some atomic lines, we have, however, found evidence of matter infall. Moreover, a narrow absorption component, which 
is particularly strong in some P-Cyg profiles, may indicate that a jet-like outflow has also been formed. 

We show that the observed emission in the [Fe II] lines and an eclipse-like event observed in November/December 2006 was probably 
caused by interactions of the expanding matter, ejected by V838 Mon in 2002, with radiation from the B3V companion. In particular, 
the observed profiles of the [Fe II] lines can be easily modelled in this scenario and allow us to estimate parameters of the system, 
such as the position of the B3V companion relative to V838 Mon and the line of sight, density in the outflowing matter, and mass lost 
in the 2002 eruption. The observed appearance of strong Ha emission, just before and during the eclipse-like event, can be interpreted 
as a result of the accretion of the outflowing matter onto the B3V companion: the accreted matter, shocked above the stellar surface, 
can be a source of extreme-UV and soft X-ray radiation capable of ionizing and exciting H in the outflow. 

Key words, stars: individual: V838 Mon - stars: late-type - stars: mass loss - stars: peculiar - stars: winds, outflows - line: profiles 



1. Introduction 

The eruption of V838 Monocerotis was discovered at the 
beginning of January 2002. As obse rved in the optical , 
. the eruption lasted about t hree months jlVIunari et all l2002at 
■ iKimeswenger et all 120021: ICrause et all l2003h . During the 
[ event, the object reached a luminosity of ~ 10^ L© . After 
developing an A-F supergiant spectrum at the maximum 
in the beginning of February 2002, the object evolved to 
[ lower effective temperatures and in April 2002 was practically 
1 unable to be detected in the optical, remaining very bright 
however in the infrared. Optical spectroscopy acquired at 
that time discovered a B3V co mpanion of the erupted object 
dMunari et alll2002bl) . lTvlendal (2005) analysed the evolution of 
V838 Mon during outburst and early decline. 

Different outburst mechanisms, including an unusual nova, 
a He-shell flash, and a stellar merger, were proposed to explain 
the eruption of V838 Mon. These m echanisms were critically 
discussed by lTylenda & Soke3 (l2006h . the authors conclude that 
the only mechanism that could satisfactorily account for the 
observational data was a collision and merger of a low-mass pre- 
main- sequence star with an ~ 8 Mq main-sequence star. 

In lKamiiiski et alj (l2009t) (hereinafter referred to as Paper I), 
a high resolution spectrum of V838 Mon acquired with the 
Keck I telescope in October 2005 was presented. In the present 
paper, we analyse and discuss the results obtained from this 
spectrum. 
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2. Evolution of V838 Mon during the decline after 
the 2002 eruption 

A few months after the eruption, discovered at the beginning 
of January 2002, V838 Mon entered a relatively calm decline 
phase. V838 Mon then resembled a very cool oxygen- 
rich (C/O < 1) supergiant slowly declining in lumino sity 
(lEvans et all l200l iTvlendaL l2005t iMunari et all l2007bh . It 
dominated the observed spectrum at green, red and infrared 
wavelengths. The blue part of the spectrum was dominated by 
the light of the B3V companion. 

Late in 2004, narrow emission lines, belonging mostly 
to [Fell], appeared m ainly in the blue part of the spectrum 
( Barsukova et all l2006). and strengthened in time jMunari et all 
I2007ah . During November-December 2006, ar i eclipse-like 
event of the B3V companion was observed (iBondl 120061: 
iMunari et all 12007 a'). At the epoch of the event, the [Fell] 
emission lines reached their maximum streng th a nd strong 
emission in B aimer lines appeared (IMunari et all 12007 al) . 

In October 2005, we obtained an optical spectrum of 
V838 Mon using the Keck/HIRES instrument. Results of these 
observations were presented in Paper I. The star V838 Mon itself 
is seen as a very cool supergiant that dominates the green and red 
parts of the spectrum. Numerous, often very deep and complex, 
molecular absorption bands are the main spectral characteristics 
of this component. All the bands are from oxides and include 
TiO, VO, AlO, ScO, and YO. The excitation temperature derived 
from the bands ranges from ~ 2500 K, identified also as the 



photospheric temperature of V838 Mon, down to ~ 500 K, 
which presumably corresponds to outflowing matter at a few 
stellar radii. The most positive (heliocentric) radial velocity 
derived from the bands is ~ 58 kms which we propose 
to be the radial velocity of V838 Mon itself. Several atomic 
lines, mostly resonance, display P-Cyg profiles. They provide 
evidence of an intense mass outflow with a typical velocity of 
~ 150 km s ' . The blue part of the spectrum is dominated by the 
spectrum of the B3V companion. The values of T^ff and log g 
derived from the spectrum agree well with those implied by the 
spectral type. The star is a rapid rotator {V sin/ ^ 250 kms"') 
and its (heliocentric) radial velocity is ~ 40 kms"'. Numerous 
emission lines were identified, mainly in the blue part of the 
spectrum. They correspond predominantly to [Fe II] and exhibit 
the same profile, which can be fitted well with a Lorentzian 
profile. The lines are centred on a heliocentric radial velocity 
of ~ 13 km s"' and have a FWHM of ~ 80 km s ' . 

iMunari et al.l ( 12007 a') propose that the appearance of the 
[Fell] emission lines and the eclipse-like event observed in 
November/December 2006 were unrelated, and that the eclipse- 
like event was caused by an eclipse in a binary system. We argue 
that both events were strongly related and produced by matter 
ejected from V838 Mon during the 2002 outb urst and reaching 
the v i cinity of the B3V companion (see also iBarsukova et al.L 
l2006HBondll2006l) . 

In principle, one can consider that because of a certain 
mechanism (e.g., dissipation of mechanical energy), the outer 
parts of the present V838 Mon wind become excited, producing 
the observed [Fell] emission lines. In this case, however, the 
radial velocity of the lines would be close to that of V838 Mon. 
This is certainlly not the case. The [Fell] lines have a radial 
velocity of 13.3 + 0.7 kms"' (Sect. 4.1.1 in Paper I), while the 
radial velocity of V838 Mon is at least 53 km s"' (see Sect.O. In 
this case, it would also be difficult to explain the observed [Fe II] 
profiles. 

All the observed characteristics of the [Fe II] emission lines 
can, however, be easily explained, if they are assumed to be 
produced by the matter ejected during the 2002 eruption that 
approaches the B3V companion. As shown in Sect. H] the 
observed line profiles as well as their radial velocities can then 
be well accounted for. The continuous strengthening of the 
lines in 2005-2006 is also easy to understand: larger and larger 
amounts of matter become excited as it approaches the source 
of excitation. Excitation by radiation from the B3V companion 
then explains why the emission line spectrum is dominated by 
Fe II. In the spectrum of a B3V star, there are enough photons 
capable of ionizing species with an ionization potential of ~ 
7 eV, e.g., Fe, Ni, while there are very few photons capable of 
ionizing the most abundant elements such as H, He or CNO. 
Finally, the idea easy explains the ~ 70 day eclipse-like event 
observed in November/December 2006 as an occupation by a 
dense cloud of the matter ejected in 2002 and now crossing the 
line of sight of the B3V companion. If this is the case, one can 
expect that the 2006 eclipse should be followed by similar events 
when other fragments of the 2002 ejecta cross the line of sight. 
Indeed, a few months later V838 Mon entered another eclipse- 
like eventQ We note that the latter would be, at least, difficult to 
explain as a phenomenon typical of an eclipsing binary system. 



' See e.g., the website of V. Goranskij 
( |http://jet ■ sao.ru/~goray/v838mon.htm^ . 




Fig. 1. A sketch of the scenario considered in modelling the 
[Fe II] line profile in Sect. E] 



3. The radial velocity of V838 Mon 

Several molecules observed in the spectrum of V838 Mon reveal 
bands formed from high excitation-energy levels (see Sect. 5 
in Paper I). They are usually indicative of a radial heliocentric 
velocity of 58 + 5 kms"' (see also Sect. 15.1b . It is reasonable to 
assume that they arise very close to the photosphere and have a 
velocity close to the stellar one. Earlier determinations of the 
stellar velocity in iKolev et all (|2002|) and iKipper et all (|2004 
implied a value of 59-62 km s"', i.e., close to our result. 

However, t h e SiO maser emission o bserve d in V838 Mon 
(D eguchi et all 120051; IClaussen et all 120071; iDeguchi et all 
20091) is at a heliocentric velocity of 71 kms"'. In the case of 
late-type stars, the SiO masers usually have a reliable measure 
of their stellar radial velocity. 

There are two possible explanations of this discrepancy. If 
the SiO maser emission really originates at the stellar velocity 
then the estimate made in Paper I may indicate that we do 
not measure molecular bands deep enough that are to reach 
hydrostatic levels of the V838 Mon atmosphere, i.e., even the 
highest excitation bands are formed in already expanding layers. 
The other possibility is that the optical spectroscopic studies 
provide a correct estimate of the stellar velocity but that the 
SiO maser is not at the stellar velocity. We note that the radial 
velocity and width of the main component of the SiO maser are 
very close to those of the CO rotational li nes observed in d iff'use 
matter in the close vicinity of V838 Mon (lKamiriskill2008h . This 
point certainly deserves further investigation. 

In the present paper, we usually adopt the radial velocity 
derived from the SiO maser, i.e., 71 km s"', as the radial velocity 
of V838 Mon. One has, however, to keep in mind that the results 
from optical spectroscopy infer a value 10-15kms"' lower. 



4. An analysis of the [Fe II] line profile 

As discussed in Sect.|2l we assume that the [Fe II] emission lines 
are formed in the matter ejected by V838 Mon during its 2002 
outburst, which is now approaching the B3V companion and 
ionized by the UV radiation of the companion. In this section, 
using simplified modelling, we show that the observed profiles 
of the [Fe II] lines can be easily explained by this scenario. 
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4. 1 . The model 

We consider a system of two stars, S i and 5 2, with a separation 
A (see Fig.[TJ. We assume that the star 5 1 is a source of a steady, 
spherically symmetric wind, that forms an expanding envelope 
with a r^^ density distribution. We allow for the possibility that 
matter in the form of clumps or filaments fills a ^ part of the 
envelope volume, i.e., ^ < 1. The wind started a certain time 
ago, so the envelope has an outer boundary, ro. The wind has 
not yet reached star S2, so ro < A. Star 52 is a source of 
ionizing radiation and emits A^ph ionizing photons isotropically 
per unit time and unit solid angle. The radiation produces an 
ionized region in the wind whose boundaries can be determined 
assuming ionization equillibrium. 

We consider ionization along a given direction from star S 2 
forming an angle a with the axis between the stars. A distance X2 
from star 5 2, to which the wind matter is ionized, can be found 
from 



A^rec = A^ph dw. 



where 

A^rec ~ 



r 



^ da> ffrec n1 dx 



(1) 



(2) 



is a number of recombinations in a cone of solid angle Au along 
the considered direction, orrec is a recombination coefficient, 
while Ni - No {A / r)^ is an ion number density (electron density 
is assumed to be equal to ion density). In the following, we 
assume that A is a unit length, so all the r and x are expressed in 
terms of A, and x is then related to r and a via 



1 + X -2x cos a. 



(3) 



In Eq. dU, xi is a distance from star 5 2 to the wind outer 
boundary along the considered direction, given by 



xi - cos a - ^cos^ a - (1 - r^). 
Equation ([TJ can be rewritten as 

_ A^ph nx\ 

Co = — T - \ -r djc, 

which, using Eq. Q, can be evaluated to give 
Co = I(x2) - I(xi), 



where 



I(x) = 



X COS 2a — cos a 
2 sin^ a 



1 



2 sin^ a 



arctan ■ 



• cos a 



(4) 



(5) 



(6) 



(7) 



Equation ^ can be solved numerically to obtain X2. Along the 
separation axis, i.e., when a - Q, xi - 1 - ro, while Eq. (|7|l 
reduces to 



/(x) = 



1 



1 



1 



l-x (l-xY 3(1 -x) 



3 ■ 



(8) 



Solving Eq. ^ for a grid of a values in the range Q < a < 
arcsinro, allows us to obtain the shape and position of the 
ionization front in the wind. Matter is ionized between the wind 
outer boundary and the ionization front. We note that the star 
separation axis is a symmetry axis of the ionized region. There 
are two free parameters in the above problem, i.e., ro and Co. 

We assume that the ionized region is isothermal, so that the 
emission line coefficient varies as A^^, and that the intrinsic line 
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Fig. 2. A model profile (full curve) compared to the observed 
profile of [Fell] 4287 A (points). The model parameters are 
discussed in the text. 



profile is Gaussian characterized by a thermal and/or turbulent 
velocity, Vi. Integrating the intrinsic line profile over the ionized 
region and taking into account the kinematic properties of the 
ionized wind, a final emission-line profile can be obtained. 
Apart from the above-mentioned parameters determining the 
ionization front, the resultant line profile depends on kinematic 
parameters of the wind and the stellar system, which are: the 
wind expansion velocity, Vwind; the velocity of star S \ (source of 
the wind) relative to the observer, Vs, and the angle between the 
stars separation axis and the line of sight, as- 

4.2. Fitting tine model to the observations 

Some of the parameters in the above problem can be estimated 
from observations. As discussed in Sect. |2l we assume that 
close to the time of the November/December 2006 eclipse, 
matter ejected during the 2002 eruption of V838 Mon reached 
the vicinity of the B3V companion. Our observations were 
completed in October 2005, so we can estimate that ro - 
0.75. During the 2002 outburst of V838 Mon, expansion 
velocities observed reached ~ 600 kms ~\ although most of 
mass loss occurred at 150^00 kms ' (|Munarietal.', l 2002al : 
ICrause et all l2003t iKipper et all l2004t Tylenda, 20051). As 



discussed in iTvlendal ( 120051) . the most intense mass loss 



occurred in March 2002, which was observed as an expanding 
photosphere of velocity ~ 270 kms"'. We therefore assume 
Vwind - 250 kms"' in the present calculations. We also assume 
that the heliocentric radial velocity of V838 Mon is = 
71 kms"' (see Sect. [3]l. Thus, there remain 3 free parameters, 
i.e., Co, as, and Vt, which can be obtained by fitting the model 
profile to the observed profile of the [Fe II] emission lines. 

A fit of this kind is presented in Fig. |2l Points show the 
observed profile of the [Fe II] 4287 A line. This is one of the 
strongest emission lines of [Fe II] in our spectrum. Its profile is 
free of blending with other lines and affected little by noise. The 
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local continuum was subtracted and the flux was normalized to 
the peak value. (Note that the profile of [Fe II] 4287 A shown in 
Fig. 12] is practically the same as the mean profile derived from 
7 [Fell] lines displayed in Fig. 4 of Paper I.) The model profile 
was obtained with the following values of the free parameters: 
Co = 035, a, = 77°, and y, = 11.5kms"'. If = 58 kms"' is 
assumed (as discussed in Sect.O, the results of the profile fitting 
remain unchanged, except as, which increases to 80°. 

We note that there is little ambiguity in determining the 
values of the above free parameters, because each of them affects 
different characteristics of the model line profile. The value of Co 
determines the size of the ionized region, so, for a fixed Vwind, 
it affects the wings of the line profile (the higher the value of 
Co, the more extended and stronger the wings). If Vwind and Vs 
are fixed, ffj defines the position of the line peak. Finally, Vi 
determines the width of the line core. 

4.3. Discussion 

The results of the above modelling of the [Fe II] profile allows 
us to estimate the parameters of the matter approaching the 
B3V companion. From the profile fitting, we have Co - 0.35. 
Assuming that the matter, expanding with velocity Vwind - 
250 kms reached the B3V companion in November or 
December of 2006, we derive a distance between V838 Mon 
and the B3V companion of A ^ 3.7 x 10'^ cm (~ 250 AU). 
Integrating the spectrum of a standard B3V star above the 
ionization potential of iron, one obtains A^ph = 2.4 x 10^^ s~' . The 
recombination coefficient of Fe I is o-rec - 2. 5 x 10~'^ cm ^ s~^ (a t 
an electron temperature of lO"* - 10'* K; see lNahar et al.L [r997h . 
From Eq. (|5]l, one then finds that No ^ 7 x 10^ ^'^'''^ cm"-', or 
A^, - 1.3x10^^"'''^ cm"^ at the outer edge of the wind envelope, 
i.e., at ro. We note that at an electron density 10'' cm"-', 
the [Fell] lines become collisionally de-excited. Assuming 
spherical symmetry, a mass-loss rate during the 2002 outburst 
can be estimated to be 

47rA2Vwmd^A?ol.4mH 

= aw^Vh ' 

where A^fc/A^h is the relative-to-H number density of elements 
with an ionization potential similar to or lower than that of Fe 
(i.e., Na, Mg, Al, Si, Ca, Ni), while mn is the H atom mass. 
Assuming solar abundances, Npe/Nu - 10""^ and Mwind - 
7 X 10^^ ^'^^gs"' ^ l^'^^Moyr"'. The most intense period 
of mass loss lasted a month (March) in 2002 (Tvlenda, 2005), 
so, according to the above estimate, it produced a ~ 0.1 ^'^^ M© 
shell. This agrees with the total mass lost in the 2002 outburst, 
which was esti mated to have a mass between 0.001 - 0.6 in 
lTvlendal(l2005l) . 

5. Mass outflow in V838 Mon 

iMunari et al. I (I2007bh noted presence of continuous mass loss 
since the outburst. They reported that the outflow observed in 
KI 7698 A has not much changed since 2002. The nature of its 
mass loss is however unclear. In this section, we analyse and 
discuss results from the profiles of molecular bands and atomic 
lines observed in October 2005 and presented in Paper I. 

5. 1 . Outflow as seen in moiecuiar bands 

As discussed in Paper I, the spectrum of V838 Mon is dominated 
by numerous, complex, and often very deep absorption lines 
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Fig. 3. The heliocentric radial velocity of individual TiO bands 
plotted against the excitation potential of the lower level of the 
band. See text for more details. 



related to molecular bands. A simple model of a stellar 
photosphere plus an outflowing homogeneous layer was used to 
identify and fit the observed structures of individual bands. In 
this procedure, it was necessary to adopt various velocities for 
the outflowing layer to reproduce different band structures. In 
general, we found that bands originating in more excited levels 
correspond to higher (more positive) radial velocities than low 
excitation bands. This effect is shown, in a more quantitative 
way, in Fig. [3] which plots the heliocentric radial velocity of a 
particular band versus the excitation energy, x, of the level from 
which the band arises. Only results for the TiO molecule, for 
which we identified the largest number of bands, are shown. 
As the excitation potential plotted in Fig. [3] we have taken 
the energy of the lower level of the most blueshifted line in 
the bandhead. Vertical error bars illustrate uncertainties in the 
velocity. Uncertainties in the energy have been taken as a range 
of energies of lower levels of lines lying within ±35 kms"' of 
the bandhead (this is the typical velocity broadening adopted in 
modelling the bands in Paper I). 

As can be seen from Fig. [3] there is a clear correlation 
between the radial velocity of a particular band and the 
excitation potential of the lower level of the band. This we 
interpret as evidence of an outward accelerated mass outflow, 
in which the observed bands are formed. The highest excitation 
bands can be formed in hotter regions, i.e., close to the 
photosphere, where the outflow velocity is low. Hence, we 
postulated in Paper I that the highest radial velocity observed 
in the high excitation bands, i.e., ~ -t-58 kms"', corresponds to 
the radial velocity of V838 Mon. The bands originating in levels 
close to the ground state can be easily formed in cold outermost 
layers outflowing with a large velocity. 

Figure |4] shows the same data as in Fig. [3] but in a 
different, more astrophysical, way. Here, instead of the excitation 
potential, we plot a parameter, logigfA) - 5040;^/rexc, where 
gf and A are the oscillator strength and wavelength of the 
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Fig. 4. The heliocentric radial velocity of individual TiO bands 
plotted against \og{gfA) - 504-0x/Te,^c- See text for more details. 

band, respectively. Within a constant factor, this parameter is 
equivalent to the logarithm of the absorption coefficient in a 
given band. As an excitation temperature, Texc, we adopted 
500 K, which is a typical value in the outflow, as inferred in 
Paper I (Sect. 5). The horizontal error bars show uncertainties 
related to the spread in the excitation energy and oscillator 
strength for the lines present within 35 km s ' of the bandhead. 

As Fig. [3] Fig. |4] can be interpreted as evidence of 
an increasing outflow velocity with the distance from the 
photosphere. Bands with a small absorption coefficient require 
large column densities to be visible in the observed spectrum. 
Thus, they must be formed deep in the outflow. More opaque 
bands are effectively formed in more distant layers and their 
profiles provide information about the kinematics well above the 
photosphere. 

From the span of the radial velocities observed in the TiO 
bands, shown in Figs. [3] and [H we can conclude that the wind in 
V838 Mon is accelerated at least to a velocity of ~ 180 kms"'. 

5.2. An analysis of the P-Cyg profiles 

The presence of lines showing P-Cyg profiles in the spectrum 
of an observed object is usually considered as strong evidence 
of an ongoing mass outflow from the object. In our spectrum 
of V838 Mon, we found about ten atomic lines exhibiting 
these characteristics (see Sect. 4.2 in Paper I). All of them but 
one originate from resonant transitions. Below we analyse the 
profiles of these lines in an attempt to investigate the physics 
and nature of the outflow. 

The terminal velocity of the outflow from V838 Mon can 
be obtained from the observed K I /17698 absorption profile (see 
Fig. 6 in Paper I). If we take Vj = 71 kms"' to be the stellar 
radial velocity (see discussion in Sect. [3]), then the terminal 
velocity of the outflow is Voo - 215 + 5 kms"'. The observed 
absorption components in other P-Cyg profiles do not reach this 
velocity, most probably because of the low column density of 



the atoms in the outer parts of the outflow. Because of the high 
oscillator strength and the high atomic abundance, the opacity in 
the K I /17698 line is at least two orders of magnitude higher than 
for any other resonant line observed. Consequently, this line can 
be observed even with relatively low column densities expected 
at high velocities. The value of terminal velocity derived from 
the KI absorption is consistent with, although slightly higher 
than, the outflow velocity estimated form the molecular bands 
of TiO analysed in Sect. 15.11 It also agrees with t he maximum 
expansion velocity derived bv lGeballe et al.l (l2007h from the CO 
bands observed half a year after our observations. 

To obtain a deeper insight into the nature of the outflow 
in V838 Mon, we performed a radiative transfer modelling 
of the P-Cyg profiles. We assume that the profiles are 
formed in a spherically symmetric expanding envelope. Our 
modelling procedur e is based on the SEI method developed by 
iLamers et alj (11987 ). and all details concerning the solution of 
radiative transfer can be found thereirQ. As the velocity field of 
the outflow, we adopt the standard y6-law, i.e.. 



w(x) = Wo -f- (1 - Wo) 1 - 




where w - V/Voo is an outflow velocity normalized to the 
terminal velocity, x = r/R-t. is a dimensionless radial distance 
from the star with a photospheric radius /?*, while wo is the 
normalized initial velocity at the stellar photosphere (taken to 
be wo=0.01 in our calculations). 

As already pointed out in Paper I (Sect. 4.2), the observed P- 
Cyg absorption components have velocity structure. A majority 
of the absorption profiles have narrow absorption components 
(NACs) that appear at different velocities. The most striking 
example of NAC is seen in the Rb I line at V/, = -82 km s"' (see 
Fig. 6 in Paper I). Some lines, e.g., those of CrI and Bal, seem to 
exhibit many NACs at once, which dominate the appearance of 
the profiles. On the other hand, there are also absorption profiles, 
e.g., those of Cal A6572 and KI /17698, which are very smooth 
and do not seem to be affected by discrete components. The 
question is whether the narrow components seen in our profiles 
are superimposed on a broad P-Cyg absorption profile formed 
in a bulk flow, or whether the entire absorption profile consists 
of several discrete and strong absorption components that are 
blended and together form a more or less continuous profile. 

As a first attempt, we assumed a distribution of optical 
thickness in the form 

T(w) = Ttot (w/wi)'^' (l - (w/wi)i/^)" , (11) 

where Tlot - ' ''"(^) is the integrated optical thickness 
and wi is the normalized velocity at which opacity vanishes 
(usually taken to be 1). The values of 7~tot, wi, ai, and 02 
were treated as free parameters in our modelling procedure. 
The radial distribution given by Eq. (fTTT i is known to reproduce 
ultraviol et (UV) r esonant lines formed in winds of hot stars 
(Castor & Lamersi 1979). Our attempts to reproduce the smooth 
profiles of Ca I, K I, and the overall profile of Rb I (neglecting 
the NAC), showed that it was impossible to reproduce 
satisfactorily any of these lines with the opacity distribution 
in the form of Eq. (fTTT i. The best set of the parameters found 



- Aware of limitations of the SEI method we have also 

i ierformed modelling in the more accurate comoving frame formalism 
Mihalas et al.Lll97a) . It confirms the results obtained with the former 
method. 
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within this approach resuhed in profiles that were too narrow 
with respect to the observed ones. 

We, therefore, modified the optical depth distribution by 
adding multiple Gaussian components, t,(w), to Eq. ( fTTT l 
parameterized by central velocities, widths, and weights 
according to which they are added. We note that, for the sake of 
convenience, the Gaussian components are parameterized in the 
velocity domain, so the corresponding distributions as functions 
of radial distance, r,(x), are in general not Gaussian (the profile 
is modified by the adopted velocity field). 

In this multi-component approach, we were able to 
successfully reproduce the observed profiles. Results for the 
Rb I /17800 line are shown in Fig. |5] At least three discrete 
components ("shells") of different widths are needed to obtain a 
profile that reproduces the observations well. These components, 
which are clearly present in all the profiles, are located at outflow 
velocities of about 50, 110, and 150 kms"' (with respect to 
V, = 71kms-i). 

We attempted to model all the prominent P-Cyg lines 
independently and, while the central velocities and widths of 
different components were found to be consistently determined 
for different lines, the relative contributions of the components 
differ considerably from line to line. This implies that the relative 
ion opacities differ considerably between consecutive "shells". 
This can be understand either as changes in excitation conditions 
or different atomic abundances (chemical composition) in the 
discrete outflow components. Although the latter seems to be 
less likely, the changing atomic fractions due to interchanges of 
atoms with molecules and/or dust grains may play an important 
role. An extreme case of this effect is observed in the resonant 
absorption lines of Ti I. Apart from the main component seen 
close to the stellar velocity (but see Sect. 15. 3b . these lines appear 
only in one of the outflowing "shells" (see Fig. 6 in Paper I). 

Our modelling allows us to place general constraints on 
the velocity field in the outflow. The adopted value of fi has 
a strong influence on the relative strengths of absorption and 
emission components in a P-Cyg profile (for a given radial 
distribution of optical thickness). We found that the outflow 
in V838 Mon can be characterized well by j6 =0.35+0.05. 
This result was obtained assuming that collisional deexcitation 
and photospheric contamination to the P-Cyg profile can be 
neglected. An influence of colUsional de-excitation on the 
emission component cannot be properly recognized because of 
unknown electron densities in the outflow. We believe that it is 
small, at least for lines with high transition probabilities (e.g., 
the lines of Rbl, KI, and Bal with an Einstein coefficient 
A„/ of the order of 10^-10^ s"'). We did not find any signs of 
photospheric contamination in the observed spectrum and we 
therefore neglected this effect in our modelling. 

The value of p derived above implies that the acceleration 
of material is very rapid and at the distance of 47?*, the 
wind has a velocity of 0.99 Vco. This steep velocity gradient 
requires very high densities of material to be able to reproduce 
strong absorption features at low velocities. Therefore, to 
reproduce the absorption component at the outflow velocity of 
50kms ' (which corresponds to r=1.03/?* for /? = 0.4), it 
was necessary to include an inner narrow shell with a very 
high optical thickness. For several lines, e.g., the KI line, it 
was even impossible to reproduce the central parts of the P- 
Cyg profile without incorporating unreasonably high values of 
optical thickness. This made us conclude that in the velocity 
field, generally well described by the yS-law, a sort of a plateau 
(or plateaus) may exist, i.e., a region of nearly constant velocity 
extended in the radial direction. These velocity plateaus, or 
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Fig. 5. Top: The P-Cygni profile of Rbl ^7800. The observed 
profile is drawn with a filled grey curve while the result of the 
SEI modelling is shown with a black line. Bottom: The mass 
loss rate corresponding to the fitted profile. The modelling was 
performed with fi - 0.4 and Vco =215 km s"^ 



non-r nonotonicities, were proposed to explain th e winds of hot 
stai-s (iPrinia & Howarthl 119861; iPuls et al.L II993 V The possible 
existence of similar features in the outflow of V838 Mon is 
a tentative conclusion, which should be examined by time- 
dependent observations of the P-Cyg profiles in the spectrum of 
V838 Mon. The derived radial distributions of optical thickness 
for different ions can be used to provide a rough estimate of the 
mass-loss rate in V838 Mon. For a gi ven stellar radius, R.^, the 
mass-loss rate can be estimated from dGroenewegen & Lamersl 
Il989h 



Mix) = t{x) w{x) x^ ^'^ 
dx 



4 nip. 



Ae flu M) ' 



(12) 



where fi is the mean molecular weight, oth is the hydrogen mass, 
Ae is the abundance of an element with respect to hydrogen, /;„ 
is the oscillator strength, Aq is the laboratory wavelength of the 
line, and the symbols nie, e, and c have their usual meanings. In 
Eq. (fT2l i. it is assumed that all the atoms of a given element are 
in the ground state and that there is no ionization nor any other 
form of ion depletion. Taking as the stellar radius = 10^ Rg 
(see Sect. 5 in Paper I), adopting /3 = 0.4, and assuming solar 
abundances, the mass-loss rate derived from the lines of Rb I 
/17800 is on average of 10"^ - lO^^Mgyr"', with individual 
mass-loss events reaching lO"-' Mq yr"' (see Fig.|5]l. 
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5.3. Signatures of an infall 

The redward broad absorption components seen in the Ti I lines 
(see Sect. 4.3 in Paper I) can be interpreted as a signature of 
an infall in V838 Mon. Similar absorption components also 
contaminate the profiles of the 5060 and 5 1 10 A lines of Fe I. 

These absorption features in the Til lines have wings 
extending out to about V/, = 120km s"'. Thus, regardless of the 
velocity of the star in the range 58-71 km s"' (see Sect. [3]), 
the red wings of these features are evidently redshifted with 
respect to the photosphere. The maximum infall velocity inferred 
from these wings is 50-60 km s"\ which agre es very well with 
the free fall velocity of an 8Mq star (Tvlenda et al.L l2005h at a 
photospheric radius of ~ 90QRq (Sect. 5 in Paper I). 

Signature s of an infa ll were previously observed in 
V838 Mon. iRushton etaTI (|2005|) observed inverse P-Cyg 
profiles of Ti I lines in their spectra obtained in the near infrared 
in December 2003. From the separation of the absorption and 
emission components, an infall velocity of ~ 20 kms ' can 
be deduced. These authors also argue that the infalling matter, 
compressed, and heated when colliding with the atmospheric 
material, provide conditions necessary to excite the SiO lines 
that they o bserve. Obser v ations of the ro-vibrational bands 
of CO by iGeballe et af] (l2007h in April 2006 revealed a 
photospheric component redshifted by 15 kms ' with respect 
to the velocity of the SiO maser. This suggests, according to the 
authors, the presence of gas infalling onto the star or contraction 
of the stellar photosphere. 

5.4. Discussion 

As mentioned above, the mass outflow in V838 Mon has been 
observed continously since the 2002 eruption. Our analysis 
confirms this statement and provides new results. The question 
that arises is what drives this outflow. It consists of cold neutral 
matter that is rich in molecules, and possibly also dust, which 
is similar to the composition of the winds of cool stars, e.g., 
objects in the asymptotic giant branch. Velocity fields of these 
cool winds are usually well described by the y6-law with p ^ 0.5, 
which is close to the value derived above for V838 Mon. These 
cool winds are believed to be dust-driven. They usually do 
not reach the high terminal velocities found for V838 Mon. 
A possible reason is that dust is destroyed in collisions, when 
the outflow is faster than about 30 kms"' (Olofsson, 2004). 
From this reason, it seems that the outflow in V838 Mon 
cannot be driven by dust only, since we observe material 
being accelerated at much higher velocities than 30 km s ' . This 
driving mechanism can, however, play a role at the base of the 
outflow, where velocities are low enough. 

One can consider a mechanism similar to that driving winds 
in hot stars, i.e., by radiation pressure absorbed in atomic 
resonant lines. However, the small number and relative weakness 
of lines of this type in our spectrum (they are typically 10 
times narrower than those in the OB stars) does not support 
the idea that these lines could be effective in driving mass loss 
from V838 Mon. More promising is the possibility that the 
outflow is driven by radiation pressure absorbed in molecular 
bands in the optical and infrared. The object is very bright in 
these wavelength ranges and the observed number and strength 
of molecular bands, as seen in our spectrum (some of them 
absorb practically all the radiation available in their wavelength 
ranges - see Paper I), clearly shows that a significant part 
of the momentum carried out in radiation is indeed absorbed 
in the outflowing matter. A process of driving winds by 



radiation pressure absorbed in molecular bands was proposed 
as a mechanism accelerating cold circumstellar envelopes (see 
J0rgensen & Johnson, 1992, and references therein). 

The analysis of the P-Cyg profiles reported in Sect. 15.21 
was completed assuming spherical symmetry and a monotonicly 
increasing outflow velocity described by Eq. ( fTOl i. Within this 
approach, the general shape of the observed profiles can be 
satisfactorily reproduced. This shows that the general pattern of 
the outflow is more or less spherically symmetric. The profiles 
however show, at least some of them, narrow structures, which 
we called NACs. The structure at Vh - -82 km s"' is particularly 
striking. In the spherically symmetric, approach the only way 
to explain the NACs is to assume that the mass-loss rate varies 
with time and produces dense shells expanding with a thinner 
wind. As can be seen in Fig. |5] to explain the strong NAC 
at Vh = -82 kms"' it is necessary to postulate a short lived 
enhancement in the mass loss rate by two orders of magnitude. 
However, if the above assumptions are relaxed, the NACs can 
be explained in other ways. One possibility is that the NACs 
are produced by mass loss being more intense in certain, more 
or less discrete, directions. For instance, the NAC at Vh = 
-82 kms"' would then be understood as a stream-like or jet- 
like outflow with a projected velocity of ~ 153 km s"' (adopting 
the radial velocity of V838 Mon of 71 km s"' ; see Sect. |3]l. Jets 
usually have counter-jets. The observed lack of any counter-jet 
in our spectrum as well as the position of the NAC close to the 
terminal velocity of the P-Cyg profiles suggests that the jet, if 
responsible for producing the NAC at -82 kms"', is not far from 
the line of sight. 

Another possibiHty is that the NAC at Vh - -82 kms"' 
has nothing to do with the presently ongoing mass outflow, but 
is produced in the matter lost during the 2002 eruption that is 
now in front of the object. This interpretation would imply the 
existence of a dense shell at a distance of ~ 1.7 x 10'^ cm 
(~ 2.5 X 10'* Ro) from the object expanding with a velocity 
of ~ 153 kms"'. However, during the 2002 eruption the matter 
was ejected w i th velo c ities ranging from ~ 100 to ~ 600 km s~ ' 
(iMunari et al.L l2002al; ICrause et all 120031; iKipper et al.L |2004 . 
Thus, the matter, or at least part of it, would have to have been 
condensed into a shell at later epochs. This can occur if matter 
ejected later catches up with slower matter ejected earlier In 
the case of V838 Mon, this mechanism does not seem to have 
worked, since the outflow velocity observed during the 2002 
eruption tended to decrease with time (Munari et al.. 2002^ 
ICrause et al.L l2003h . A NAC at Vh = -90 km s"' is seen in 
our spectrum in the Til 5147 and 5152 A lines, which do not 
display P-Cyg profiles (see panel c in Fig. 6 in Paper I). The 
close radial velocity and line width suggest that this component 
is produced in the same enviroment as the NAC in the P-Cyg 
profiles. If this is the case, our interpretation presented above 
would have to be excluded. The Ti I lines arise from the levels 
lying 2.3-2.4 eV above the ground state. Thus, they must arise 
in a rather warm matter (of temperature comparable to that in the 
photosphere), while the matter lost in 2002 and now at a distance 
of ~ 30 photospheric radii is expected to be cold (~ 500 K, if 
estimated from pure geometric dilution of the radiative energy 
density). In conclusion, we consider that the interpretation that 
the NAC at Vh - -82 kms"' is related to matter lost in 2002 is 
less probable than the other two, i.e., varying mass-loss rate in 
the present wind or jet-like ouflow. 

Spectroscopic data allowing analyses of the time-dependent 
behaviour of the P-Cyg profiles would help to discriminate 
between interpretations involving varying mass-loss rate, jet- 
like outflow, or distant matter lost in 2002. In the first case. 
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the NACs are expected to migrate in the profile with time, 
while in the two other cases the NACs should be rath er stable 
features. It is worth noting that iGeballe et alj ( |2007[) in their 
high-resolution spectra obtained in the K band in April 2006, 
i.e., half a year after our spectroscopy, indentified several 
velocity components in the CO absorption line profiles. Apart 
from the photospheric component (mentioned in Sect |5.3l l the 
authors found components at expansion velocities of 15, 85, and 
150 kms~'. The last feature is practically at the same velocity 
as our NAC at Vh = -82 kms"' suggesting that this is a 
persistent component. However, the feature in the CO lines is 
much wider than the NAC, so the conclusion is unclear. The two 
other components in the CO lines do not match other NACs in 
our P-Cyg profiles, which are seen at the expansion velocities of 
50 and 110 km s"' (see Sect.l7!2]). This suggests that these NACs 
are transient and produced by varying mass-loss rates. 

6. Summary and discussion 

Our analysis of both the positions of the molecular bands 
(Sect. 15.11 ) and the P-Cyg profiles (Sect. I5.2l i have shown that 
V838 Mon loses matter with a terminal velocity of ~ 215 km s"' 
at a rate of 10"^ - 10"^ M0yr"'. Thus, the object loses of the 
order of 10 in the form of the kinetic energy of the wind. 
This is small compared to the radiation luminosity of the object, 
which is ~ 2.4 x 10"* Lq (see Sect. 5 in Paper I). However, the 
momentum carried out in the wind is comparable to that in the 
radiation. Thus, the wind in V838 Mon is likely to be radiation 
driven, probably due to absorption in the molecular bands, as 
discussed in Sect. 15.41 

The kinematic pattern of the matter in the atmosphere of 
V838 Mon is very complex. Both expansion and wind outflow 
dominate. Most of the outflow is more or less spherically 
symmetric as suggested by our modelling of the observed P-Cyg 
profiles (Sect. |572] i. However, as discussed in Sect. 15 .4[ it is quite 
possible that we also observe a jet-like outflow. As shown in 
Sect. 15.31 certain regions, perhaps initially accelerated but not 
sufficiently to allow them to leave the object, also fall back to 
the photosphere. We note that simultaneous signatures of infall 
and outflow have been observed in the winds of hot stars (e.g., 
iHowk et aLl l2000l) and often in protostars. A possible scenario 
for creating coexisting mass outfl ow and infal l in asymptotic 
giant branch stars was discussed in Soker (2008) and may also 
work in the remnan t of the V838 M on outburst. 

As discussed in lTvlendal (l2005l) . the evolution of V838 Mon 
in the post-outburst state is probably dominated by gravitational 
contraction, i .e., is closely rel a ted to that of protostars. If, 
as argued in iTvlenda & Sokeil (2006), the 2002 outburst of 
V838 Mon was produced by a stellar merger, then a significant 
amount of angular momentum sho uld have been stored i n 
the remnant of the event (see also ISoker & Tvlendal |2007|) . 
Gravitational contraction may then result in the flattening of the 
envelope, eventually leading to the formation of an accretion 
disc (as in protostellar objects). The possibility discussed in 
Sect. 15.41 namely that we observe a jet in the spectrum, 
if confirmed, may imply that a disc has already formed in 
V838 Mon. 

As shown in Sect. |4] the observed profile of the [Fe II] lines 
can be easily understood as having been produced by matter 
lost by V838 Mon in the 2002 eruption and now excited by 
the radiation of the B3V companion. As the matter approaches 
the companion, a larger and larger proportion of the flux from 
the star can interact with the matter causing an increase in 
the [Fell] emission, as observed. Indeed, the flux in the 20 



strongest [Fe II] lines in our spectrum obtained in October 2005 
in creased by a factor of ~ 3.5 in the measurements completed 
by iMunari et al. (2007a) in December 2006. At the moment of 
our observations, the outer boundary of the expanding envelope 
was presumably at a distance of ~ 0.75 times the separation 
between V838 Mon and its B3V companion. Thus, only ~ 17% 
of the radiation flux of the B3V companion should then interact 
with the matter At the moment of the eclipse-like event in 
November/December 2006, the matter presumably reached the 
closest vicinity of the companion and the portion of the stellar 
radiation capable of exciting the matter in creased to 50 - 100% . 

Apart from the [Fe II] emission lines, Munari et al. I (I2007ah 
also observed strong Balmer lines in emission. In particular, they 
discussed in detail the profile of Hot observed in December 2006 
- April 2007. No Ha was seen in our spectrum in October 2005. 
However, in the absorption cores of higher members of the 
Balmer series (in the spectrum of the B3V companion), we have 
detected weak emission components (see Sect. 4.1.2 in Paper I). 
In particular, the emission component in the Hfi line was 
estimated to have an observed flux of ~ 3 x 10"'^ erg cm"^ s"' . 
We note that the profile of the H/3 emission component (see 
Fig. 5 in Pap er I) was rem arkably similar to that of the Ha line 
observed bv , Munari et al.l (l2007a ) in December 2006. Moreover, 
the parameters of the emission component in Hj6 in our spectrum 
and of the Ha emission component in Munari et alj (l2007al) 
(radial velocity and FWHM) are close to those of the [Fell] 
lines derived in Paper I. This implies that all these lines probably 
originate in the same region, i.e., in the part of the V838 Mon 
envelope (ejected in 2002) approaching the B3V companion. 
The point is, however, that while in the spectrum of a B3V 
star there are enough photons capable of ionizing iron and 
exciting the [Fe II] lines, there are very few photons, capable of 
ionizing hydrogen. Between October 2005 and December 2006, 
the flux in the Wfi emission also increased by a factor of ~ 7, 
while that in Ha increased by at least two orders of magnitude. 
This large increase cannot be explained by a purely geometric 
effect, as discussed above in the case of the [Fell] lines. 
Therefore, an additional source of excitation of the Balmer lines 
is necessary, a source, whose intensity significantly increased 
between October 2006 and December 2007. We suggest that this 
is provided by matter accretion from the outflowing V838 Mon 
envelope on the B3V companion. We now complete simplified 
estimates for this scenario. 

We assume that matter flows with a velocity, Vwind, in the 
vicinity of a star of mass, Mb3v- An accretion radius, race, can be 
estimated from the standard formula (see e.g.. lFrank et al.Ll2002l) 

r^^^ . . 3 X lO'^cm ^ A?>R^, (13) 

^wind 

where we adopted Mb3v - 7Mq and Vwind = 250 km s"'. An 
accretion rate can then be estimated from 

2 

- 7?%TA^wmd - 1 X lO'V'^'gs"', (14) 

(471 A^) 

where we used the values of M^md ™d A estimated in Sect. 14.31 
This implies an accretion luminosity, 

Lace = 1 - (15) 

"B3V 

where /?b3v = A-.SRq has been assumed. This luminosity is 
small compare d to the intrinsic lumin osity of a standard B3V star 
(~ 1.9x10^ Ln. lSchmidt-Kaleitll982l) . Thus, no significant eff^ect 
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is expected to be observed in the global energetics of the B3V 
companion, e.g., in the visual brightness. However, after being 
dissipated via an accretion shoclfland thus producing a hot layer 
above the star surface, the accretion luminosity is expected to be 
emitted mainly in the far-UV and/or soft X-rays. This can affect 
the ionization of the expanding envelope of V838 Mon, i.e., 
produce significant ionization of hydrogen. The accretion from 
the expanding matter is expected to become effective just before 
and during the November/December 2006 eclipse. We suggest 
that this is the reason for the rapid increase in emission of the 
Balmer lines, particularly in Ho-, observed in the second half of 
2006. We note that the luminosity in the Balmer emission lines, 
dominated by Ha, as observed bv i Munari et al.l (|2007a) during 
the eclipse, was ~ 3.5 L© (assuming Eb^v - 0.9 and a 8 kpc 
distance to V838 Mon). Thus, the accretion, as estimated above, 
can provide enough energy to explain the observed luminosity 
in the Balmer lines. 
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